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ABSTRACT 

We study the average X-ray and soft y-ray spectrum of Cyg X-1 in the hard spectral state, 
CN ■ using data from INTEGRAL. We compare these results with those from CGRO, and find a 

good agreement. Confirming previous studies, we find the presence of a high-energy MeV tail 
beyond a thermal-Comptonization spectrum; however, the tail is much softer and weaker than 
[tJ ' that recently published by Laurent et al. In spite of this difference, the observed high-energy 

l>H ' tail could still be due to the synchrotron emission of the jet of Cyg X- 1 , as claimed by Laurent 

•-H ; et al. 

f~| ' In order to test this possibility, we study optically-thin synchrotron and self-Compton 

0^', emission from partially self-absorbed jets. We develop formalisms for calculating both emis- 

JL ■ sion of the jet base (which we define here as the region where the jet starts its emission) and 

^ ' emission of the entire jet. We require the emission to match that observed at the turnover 

"^ , energy. The optically thin emission is dominated by that from the jet base, and it has to be- 

j^ ■ come self-absorbed within it at the turnover frequency. We find this implies the magnetic field 

strength at the jet base of Bq oc Zq, where zo is the distance of the base from the black-hole 
centre. The value of Bq is then constrained from below by the condition that the self-Compton 
emission is below an upper limit in the GeV range, and from above by the condition that the 
Poynting flux does not exceed the jet kinetic power This yields Bq of the order of ~ lO'* G 
^o ' and the location of the jet base at ~ 10^ gravitational radii. Using our formalism, we find the 

tH- ' MeV tail can be due to jet synchrotron emission, but this requires the electron acceleration at 

(SJ , a rather hard power-law index, p ^ 1.3-1.6. For acceleration indices of p > 2, the amplitude 

^J • of the synchrotron component is much below that of MeV tail, and its origin is likely to be 

^—^ I due to hybrid Comptonization in the accretion flow. 

^^ I Key words: accretion, accretion discs - radio continuum: stars - stars; individual: Cyg X-1 

- stars: individual: HDE 226868 - X-rays: binaries - X-rays: stars. 

>: 

Cu ■ 1 INTRODUCTION magnitude as well as much softer, with a ^ 2. The weakness of 

^ ,^ , . , . , , „ ,-,,,,,,,■ T the tail is also confirmed by an independent study of the average 

Cyg X-1 IS an archetypical and well-studied black-hole bina^'. In j^jeGRAL SPI spectrum ^urdain. Rogues & Malzacll20i2h . and 



> 

0^ 



it is in agreement with the CGRO results of M02. 



the hard spectral state, its hard X-ray /soft y-ray spectrum has been 

known to be hard below the maximum [in the EF(E) representa- 

tion] at -100-200 keV. Above it, the spectrum shows a sharp cut- ^^ '"i"*" ^^^°'^' *e high-energy tail can be well fitted by 

ff u- u .1, ■ f 11 J u u- 1, I 1 J . emission of non-thermal electrons forming a high-energy tail be- 

ott, which then is rollowed by a h igh-energy tail, measured up to a & t. bj 

few MeV JMcConnell et aljB . her eafter M02). y™'^ '^^ Maxwellian distribution, responsible for the bulk of the 

We consider here the finding of Laurent et al. I JJoTl (here- ^-'""y ^^^^^^°^- On the other hand, the overall spectrum can also 

after LI 1) of the tail at 0.4-2 MeV being very strong and hard, *"= fi"«'^ ^y P'^'^^y ^'^^^l Comptomzation (presumably m a hot 

■ ,, ,, .ij c n^.noTu.i J accretion flow) and emission of power-law electrons with an expo- 

with the energy spectral index ot a == 0.6 ± 0.2. The tail measured / i- f 

uTii c jiu» 1 i-j u- u ij- i» nential cutoff", which electrons could be located within the jet. 

by Lll was round to be strongly polarized, which would imply it ■' 

is due to synchrotron emission of the jet in this system. We calcu- I" order to test the jet origin of the tail, we study emission of 

late the average hard-state INTEGRAL spectra from the ISGRI and partially self-absorbed jets, such as that present in the hard state of 

PICsIT detectors, and although we do find a high-energy tail above Cyg X-1. We develop two formalisms. In one, a one-zone model 

about the same energy as Lll, it is weaker by about an order of of *e jet base is developed. In the other, we model the emission 

of the entire jet (reformulating the model of iBlandford & Konigil 

I1979I hereafter BK79). We find that the optically-thin emission is 

* E-mail: aaz@camls;.edu.pl dominated by the jet base. In both approaches, we relate the mag- 
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Figure 1. The hard-state radio to y-ray spectrum of Cyg X-1. The five radio/mm fluxes are from lFender et alj J2OOOI) . The two IR broken power-law spectra 
show the jet component for o bservations 1 and 2 of Rl 1. The dott ed line shows an d = radio spectrum extrapolated up to 0. 1 eV. The cyan circles and error 
bars show the total IR fluxes jPersi et al.ll980l : lMirabel et al.ll996) . which are mostly from emission of the companion and its wind. The black X-ray and soft 
y-ray data points above 20 keV show the IBIS and COMPTEL spectra. The spectrum below 20 keV, from BeppoSAX, represents a typical hard state spectmm 
(which is absorbed by an intervening medium). The dot-dashed magenta curve shows a fitted model of the unabsorbed spectrum. The y-ray upper limit is from 
AGILE iSabatini et allioiOi) . See Sections[2jg]for details. 



netic field at the jet base to its height along the jet. We apply our 
model to Cyg X-1, using available measurements and upper limits. 
We find that the observed turnover energy, the flux at that energy, 
the MeV tail measurements, an upper limit on the emission in the 
GeV region, and a measurement of the kinetic jet power allow a 
relatively precise determination of both the height of the jet base 
and its magnetic field. 



2 NOTATION AND THE PARAMETERS OF CYG X-1 

Hereafter, we use the following dimensionless parameters, 



£■' 



^•h 



Mc^ 






(1) 



where E' is the photon energy in the jet comoving frame, z is the 
height along a jet (from the centre of the black hole in the observer's 
frame), R^ = GM/c^ is the gravitational radius, G is the gravita- 
tional constant, M is the black hole mass, M is the mass accretion 
rate, Pj is the sum of the total kinetic powers of the jet and coun- 
terjet, Lg = 4nGMmpC/crj is the Eddington luminosity, ctj is the 
Thomson cross section, and nic, nip and are the electron and pro- 
ton mass, respectively. We use then the symbol E for the observed 

photon energy. 

We adopt the best-fit values of lOrosz et"ai] ( 1201 ih of M ^ 
15M0 (where M© is the solar mass) and the inclination of the 
normal to the binary plane with respect to the line of sight of 
; = 27°. We assume that the jet and the disc have the same in- 
clination. \^_use the best-fit value of the distance to Cyg X-1 of 
iReid et al.l j201lh of D = 1.86 kpc. The separation between the 



star s is a - 3 X 10'^ cm, and the stellar luminosity i s L, ^ 10^^' erg 
s' JCaballero-Nieves et al.ll2009l:IOrosz et"ani201 lb . 

The opening angle of the jet in Cyg X-1 on the length 
scale of ~ 10'"' cm has been co nstrained by VLBA and 

VLA observations to 0j ^ 2^ ('Stirling et al.' '2OOI I), and 

its ve locity, to y6| ^ 0.6 1 Stirling e t al. 2001; Gleissn er et al.l 
l2004l ; lMalzac. Belmont & Fabianll 20091) . which corresponds to the 
Lorentz factor of Fj = (1 - ff'T^^^ - 1.25. The kinetic jet power 
in the hard state of Cyg X-1 has been estimated as P^ =^ (1- 
3) X 10 ^^ erg s~' from the optical nebula presumably powered by 
the jet dGallo et al. 2005; Russell et al. 2007). For M = 15Mq, 
Le - 2 X 10^' erg s"', which gives pj ~ 10"^. On the other hand, 
ra ^ 0. 1 is obtained for the observed average X-ra y luminosity of 
Lx - 10'^ erg s ' in the hard state of Cyg X-1 dZdziarski et al.l 
I2OO2I) using an accretion efficiency of 0.05. 

The pos i tion a ngle of t he radio jet on mas scale found by 
IStirling et aO ( l200lh (see also iRushton et alj|201 ll) is -(17°-24°). 
However, the sign o f these values was misprinted as positive in 
IStirling et al.l ( I2OOII) . The actual position angle, conventionally 
measured from the n orth to the east, is negative, as stated, e.g., in 
ISpenceretaLl ( l200ll) . 

A broad-band spectrum of Cyg X-1 in the hard state is shown 
in Fig. [T] In the radio/mm range, we s how the fluxes at 2 .3, 8.3, 
15, 146 and 220 GHz (from Table 1 of iFender et alll2000l) . These 
fluxes can be joined by an cc = spectrum, which extends up to 
the turnover energy at E, == 0. 1 eV in the IR. The two 0.05-0.24 
eV broken power-law spectra show the fitted jet component for ob- 
servations 1 and 2 (case 2, model c) of iRahoui et alj ( l201lh (here- 
after Rll), with the fitted values of £, = 0.12 eV and 0.06 eV, 
respectively. The cyan crosses and error bars show the measure- 
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ments of IPersi et af] ( Il980l) and iMirabel et al] ( Il996l) of the total 
IR emission, which is dominated by the companion and its wind. 
Note that the fitted jet component represents a relatively small frac- 
tion of the total IR emission, especially at its high-energy end. 
The black 20 keV-5 MeV data points show the average ISGRI 
and PICsIT spectra obtained by us (see below), and the average 
0.7-5 MeV COMPTEL spectrum from M02. The spectrum shown 
at 0.5-20 keV is an e xample of a typical hard-state spectrum, ob- 
served by BeppoSAX JDi Salvo et al.l200lh . Note that this spectrum 
is strongly absorbed at low energies by an intervening medium; the 
magenta dot-dashed curve shows a model fit of the intrinsic spec- 
trum. The 0.1-3 GeV upper limit is from AGILE JSabatini et al.l 
|2010|) . with the integral photon flux of < 3 x 10"^ cm"^ s"' (2cr). 
Given the large width of the energy band with a strong dependence 
of an implied monochromatic upper limit on the s pectral shape, 
we use here the same EF(E) limit as that of fig. 3 of lSabatini et al.l 
1 120101) . This approximately corresponds to the monochromatic flux 
of £ Fy = 6 X 10"' cm"^ s"' in the middle of the observed range 
at Ey = 0.5 GeV. The present upper limit is a factor of a few belo w 
that of the CGROfEGRET (given in lZdziarski & Gierlinskill2004l) . 



3 THE X-RAY AND SOFT y-RAY SPECTRA 

3.1 Analysis of the INTEGRAL data 

We study INTEGRAL data from the IBIS instrument. We consider 
separately its two detectors, the ISGRI and PICsIT. We use data 
from the INTEGRAL revolutions from 22 to 877 (MJD 52626- 
55186, the years 2003-2009), which are available in the INTE- 
GRAL Science Data Centre (ISDC). All public data with the oflF- 
axis angle < 15° were selected (a < 9° condition gave virtually the 
same spectra). The ISGRI and PICsIT data are used in the 19-500 
keV and 0.25-5.4 MeV ranges, respectively. 

We obtain both the average spectrum from all the above 
data, and the average spectrum in the hard spectral state. The 
former is still dominated by the hard spectral state, which is 
due to Cyg X-1 having been mostly in it since the launch of 
INTEGRAL until 2010 June. Thus, the diiferences between the 
two sets are relatively minor As the hard-state intervals, we 
take MJD 50350-50590, 50660-50995, 51025-51400, 51640- 
51840, 51960-52100, 52565-52770, 52880-52975, 53115-53174, 
53540-53800, 53880-55375, following Zdziarski (2012). In addi- 
tion, we also use the data from MJD 53188-53215, during which 
the spectra appeared typical to the hard state. 

Our data analysis procedure is as follows. The ISGRI data 
have been reduced using the Offline Scientific Analysis (osa) 9.0 
provided by the ISDC (C ourvoisieret al. 2003) . with the pipeline 
parameters set to the default values. For the spectral extraction we 
took into account also 22 strongest ISGRI sources in the Cyg X- 
1 neighbourhood. The PICsIT spectra were prepared using a non- 
standard software, optimized for handling Poisson-distributed data 
affected by a strong and variable background ( Lubiriski..2009.) . We 
have added a 1 per cent systematic error to each of the ISGRI and 
PICsIT data sets. 



3.2 Average spectra 

We have fitted our IBIS spectra for both all the observations and 
those in the hard state only with the model of hybrid Comptoniza- 
tion, EQPAiR, obtaining results similar to those of M02. The de- 
convolved spectra for the two data sets are shown in Figs. |2la) 



and (b), respectively. Since the fit results are similar for both data 
sets, and since the former contains also some contribution from 
soft/intermediate states, we describe the fit results only for the latter 
below. We note that the overall flux calibration for the PICsIT dif- 
fers slightly from that for the ISGRI, yielding the PICsIT flux about 
1 .3 times higher than that of ISGRI in the overlapping energy range. 
In Fig.|2ta), the PICsIT spectrum is shown at its original normaliza- 
tion, whereas in Figs.|2jb)-[3]it is plotted renormalized to the level 
of ISGRI. We note that such differences between the noiTnalization 
of different instruments are common, e.g., a similar difference was 
present for the PCA and HEXTE detectors onboard RXTE. 

Fig.|2ta) compares then our IBIS spectrum with that published 
in LI 1. The latter spectrum is from the ISGRI up to ~400 keV, and 
at higher energies from the Compton mode, which utilizes both IS- 
GRI and PICsIT. We see that those fluxes are substantially above 
ours, by about 30 per cent at 100 keV, and up to a factor of several 
above I MeV. However, it turns out (Ph. Laurent, private communi- 
cation) that the published spectrum does not represent the average 
Cyg X-1 spectrum from the IBIS, and it does not correspond to the 
data used to calculate the polarization. The corrected average spec- 
trum (Ph. Laurent, private communication) is shown in magenta 
and cyan in Fig. |2ja). Its ISGRI part agrees with our ISGRI spec- 
trum almost perfectly. The Compton-mode spectrum still shows a 
disagreement with our PICsIT spectrum at ^ I MeV, but now it is 
at a level of ~(2-3)cr, compared to several cr discrepancy for the 
published spectrum. 

In Fig.|2la), we als o show t he average spectrum from the IN- 
TEGRAL SPI detector of ljourdain et al.. C201Z) . This spectrum was 
obtained using all usable SPI Cyg X-1 data, which time interval 
spans 52797-55184, and which, similarly to our average, also in- 
cludes some fraction of soft/intermediate-state data. Although there 
is an overall agreement with our results, the SPI spectrum is above 
ours at £■ ^ 30 keV, and shows a much harder decline up to 700 
keV. This may be partly due to the ISGRI calibration with re- 
spect to the Crab spectrum, assuming that spectrum above 100 keV 
has a = 1.35, whereas c urrent estimates give a = 1.24 ± 0.02 
jjourdain & Roauesll2009l) . This also seems to explain our PICsIT 
spectrum being harder than that from ISGRI at S 150 keV. We also 
see that the 2cr upper limit from SPI at 1.1-2.2 MeV is above our 
PICsIT flux but still somewhat below the corrected Compton-mode 
flux (of Ph. Laurent, private communication). 

Fig. Ob) compares then our INTEGRAL/lBlS average spectra 
for the hard state with the averaged ones for 1991-1994 measured 
by CGRO (M02). Given uncertainties in the instrument calibra- 
tion, we have not applied any correcting factors to the BATSE and 
COMPTEL data, which were applied in M02 in order to normalize 
those data to the OSSE spectrum. We see that our IBIS spectrum is 
very similar to that from BATSE except for a slight difference in the 
normalization. At E > 0.4 MeV, the PICsIT, COMPTEL, BATSE 
and OSSE spectra all agree well. 

The EQPAIR model is described in detail in ICoppil ( 1 19991) and 
iGierliiiski et al.l ( 119991) . Our best-fit parameters for the hard state 
are: the ratio of the power supplied to the Comptonizing plasma 
to that in the seed photons irradiating it of 6.8 ± 0.2, the fraction 
of the power supplied to the non-thermal electrons in the plasma 
of 0.16 + 0.01, the Thomson optical depth of ionization electrons 
of T; = 0. 79 ± 0.07, the relative strengt h of reflection (using the 
method of 'Magdziarz & Zdziarskil|l99 5') of n/ln = 0.25 ± 0.03. 
There is a relatively strong e* pair production in this model, with 
the total Thomson optical depth including pairs of 1 .44, and at the 
equilibrium temperature of e* of kT^ - 69 keV (note that these two 
quantities are not parameters of the model but instead result from 
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Figure 2. Comparison of average 20 keV-5 MeV spectra of Cyg X- 1 . (a) 
The spectra averaged over all INTEGRAL observations from 2003-2009. 
Tlie black and blue spectra are from our analysis of the data from the ISGRI 
and PICslT, respectively, fitted together by hybrid Comptonization. The red 
spectra are from the ISGRI and the Compton mode (ISGRI+PICslT), as 
published by Lll. The magenta and cyan crosses are the coiTected ISGRI 
and Compton-mode spectra, respectively (Ph. Laurent, private communica- 
tion). Each of those two pairs of spectra have been fitted by thermal Comp- 
tonization and a power law. The green spectr a are from the SPI, fitted by 
two thermal-Compton components (Jourdain etal .1120121) . (b) The average 
spectra in the hard spectral state from INTEGRAL and CGRO. The black 
and blue spectra are from the ISGRI and PlCslT (with the PICsIT data 
renormalized to the level of the ISGRI). The magenta, cyan and red sym- 
bols show the CGRO OSSE, BATSE and COMPTEL spectra, respectively, 
fitted together by hybrid Comptonization (M02). See Section[3]for details. 
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Figure 3. The average 20 keV-5 MeV spectra of Cyg X- 1 from our analy- 
sis of the ISGRI (black) and PICslT (blue) in the hard spectral state fitted 
by (a) hybrid Comptonization and reflection and (b) thermal Comptoniza- 
tion, reflection, and an e-folded power law. For comparison, the COMPTEL 
spectrum from M02 (not taken into account in the fits) is shown in red. The 
short and long dashes and the solid curves show the thennal Compton scat- 
tering, Compton reflection, and the total model, respectively. The dotted 
curve in (a) shows the Compton emission from the non-thermal e* plus the 
annihilation spectrum from the e* pairs, and in (b), the e-folded power-law 
component. See Section[3]for details. 
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the fitted and assumed parameters). We have assumed the compact- 
ness parameter of 10 of the blackbody seed photons at the temper- 
ature of 0.2 keV. The amount of pair production depends on the 
former , and the dependence o n the latter is rather weak, see discus- 
sion in lGierlinski et al.l 1 119991) . The model and its components, i.e., 
scattering by thermal electrons, scattering by non-thermal electrons 
plus pair annihilation (which is due to pairs produced by the non- 
thermal part of the spectrum), and Compton reflection, are shown 
in Fig.[3ja). An annihilation hump is seen in the second component. 
The decomposition of the Compton spectrum into the thermal and 
n on-thermal components ha s been done using the method described 
in lHannikainen et alj 120051) . We have assumed that the model has 
a 1 per cent systematic error. We have obtained a relatively good 
fit, ;('^/dof = 94/58. The model is shown by the solid curve in Fig. 
Ob). Note that although we have not fitted the COMPTEL data, our 
model agrees well with them. 

We have then tested a two-component model in which the 
main part of the spectrum was fitted by thermal Comptonization 
without any non-theimal tail, and the observed tail was attributed 
to a power law with an exponential cutoff. We have found that the 
energy index, a, of the power law could not be constrained as the 
same^' was found in a wide range of a. However, since a possible 
physical inteipretation of this component is a synchrotron contri- 
bution from the jet, we have fixed it at o- = 0.6, for which this 
power law extrapolated down to the turnover energy of £t - 0. 1 
eV connects to the extrapolation of the optically thick jet emission 
with a ^ 0, see Fig. \T\ We have obtained the e-folding energy of 
E^ = 1.2±0.3 MeV and the 1-keV normalization of 0.11 cm"^ s ', 
at^^/dof = 96/57, i.e., only slightly worse than that of the hybrid- 
Compton model. The thermal plasma parameters are: the power 
ratio of 6.5 ± 0.2, Tj = 1.52 ± 0.03 (with virtually no pair produc- 
tion, and with kT^ = 69 keV), and Q/ln = 0.26 + 0.03. We see 
that the two-component model predicts a hump at ~2 MeV some- 
what more pronounced than the hybrid-Compton model. The two 
components intersect at 410 keV Thus, if the e-folded power law 
were strongly polarized, it could explain the polarization measured 
above 400 keV by LI 1. 



4 MODELS OF JET EMISSION 

We study both a one-zone model of the jet base region and a jet 
model of the entire jet, and apply it to Cyg X-1. For our purposes, 
we define the jet base as the location at which the jet emission starts 
(assuming continuous jet emission, as in the model of BK79). This 
is likely to differ from the actual base where the jet is formed but 
the dissipation has not yet started. 

As we confirm in Section 14.21 the base provides most of the 
optically-thin jet emission. The jet base becomes optically thin 
at the turnover energy, Et. The value of this energy and the syn- 
chrotron flux at Et in Cyg X-1 have been measured by Rll. This 
allows us to determine the normalization of the distribution of the 
emitting electrons as a function of the magnetic field. The syn- 
chrotron emission has to be self-absorbed within the jet base at 
Et, which imposes a relationship between the magnetic field, Bg, 
and the location of the base, zo- Then the magnetic field strength is 
constrained from below by the requirement that the self-Compton 
emission in the y-ray energy range does not exceed the observa- 
tional upper limit. On the other hand, the magnetic energy flux can- 
not be greater than the observed jet kinetic energy. 



4.1 One-zone model of the jet base 

The radio up to far IR emission of Cyg X- 1 is flat in the dF/dE 
representation, oc E", with a - (see Fig. [TJ. This is usually 
interpre ted as a partially synchro tron self-absorbed jet emission 
(BK79; iFalcke & Biermanrj 1 1995 1). The emission at a given en- 
ergy is self-absorbed from the base of the jet up to the height of 
z o^ V"', and it is optically thin at higher z- The jet becomes op- 
tically thin at all z above £,. For a magnetic field of 6 oc z"', 
conserving the magnetic energy flux (or for steeper dependencies), 
the optically-thin emission ai E > Et is dominated by the emis- 
sion from the base of the jet. This follows, e.g., from the syn- 
chrotron flux, E(E) oc g(p+i)/2^ where p is the index of the elec- 
tron power law. For B oc z^\ F{E) oc ^-(p+1)/2 jj^g number of 
emitting electrons for a constant- velocity jet per In z is oc z- Thus, 
dF(£)/dlnz oc z^", where a = (p - l)/2 is now the index of the 
optically thin synchrotron radiation. Integrating over z, we see that 
the first logarithmic interval of the jet length measured from its base 
provides about a half of the total emission for a = 0.5. Thus, in this 
Section we adopt a one-zone approximation, with the emission re- 
gion at z + z/2 (in the observer's frame) along the jet. The validity 
of this approximation for optically-thin emission is confirmed by 
our calculations of the synchrotron spectrum from the entire jet in 
Section|42l 

We consider a broken power-law steady-state electron distri- 
bution per unit volume. 



A.(.).^exp(^)x{;;V. 



1 ~ To < r ^ Tb; 



(2) 



where K is the normalization, y is the Lorentz factor, and the cool- 
ing break, jf,, and the cutoff, 7,^ (assumed » yj,) are discussed and 
determined in Section 1431 below. In Sections 14.1114.21 we assume 
that the distribution of electrons responsible for the emission in a 
region around the turnover energy corresponds to y < yb, where we 
can also neglect the high-energy cutoff, i.e., N(y) = Ky^''. 

The synchrotron emission coefficient (i. e., energy production 
rate p er unit volume and per unit solid angle, iRvbicki & LightmanI 
1 19791) from isotropic relativistic electrons in a unit volume at the jet 
frame for a gen eral form of A^(y ) can be approximately described 
as (see, e.g. Zdziarski et al.ll2012l hereafter Z12), 



js(e) i 



1/2,1/2 



edhs CTjcBUB/B^tY'^e 



dedQ 



48;r2 



-N 



B 



(3) 



where B^r = 2nmlc^/eh is the critical magnetic field, h is the 
Planck constant, e is the electron charge, and h denotes a pho- 
ton production rate per unit volume. The formulation of the syn- 
chrotron coefficients in terms of ctj and fi^r expresses the close cor- 
responde nce between the synch r otron and Compton processes, dis- 
cussed in lBlumenthal & Gould ( Il970l) . The angle-integrated emis- 
sivity per unit volume is 4njc, and the power emitted in all direc- 
tions by a volume V in the jet frame is AnVj^. (We note that the 
similar formalism used in Z12 employed rates integrated over the 
source volume.) For power-law electrons, we have approximately 



Js(e) ^ 



CiOTTcKBi I B 



48;r2 



Be 



£ - , 



(4) 



where Ci = 1 in the delta-function approximation of equa- 
tion ([3}, w hereas averaging over the pitch angle gives (cf. 
Ijones.O'D ell & Stein 1974) 



C, 






(5) 



© 2012 RAS, MNRAS OQO.fnTBl 



A. A. Zdziarski et al. 



where T is the gamma function, and Ci = 1 for p = 3. 

The synchrotron self-absorption coefficient averaged over the 
pitch angle for a power-law electron distribut ion can be expressed 
as ( 



tch angle for a power-law electron distribut io 
; (cf. lJones et al.lll974lzdziarski et al.ll2012l) . 



CjnojKl B 



lar 



lar Be, 



B 



C2 = 



3^r(^)r(^)r(^) 



, (6) 



(7) 



where Uf is the fine-structure constant and C2 - 0.9996 for p = 3. 
The second equality in equation {6} gives a monochromatic approx- 
imation for as, which is almost accurate f or p = 3. 

The self-Compton emission is approximately given by (e.g., 
Z12), 



jsc(e) ■ 



1 1 /' 



Sn 



Jo 



TT. A'L - deo, (8) 



where ns + «sc is the density of the synchrotron and self-Compton 
photons. 



(«s + nsc)(fo) - { [jsi^o) + Jsc(eo)J ■ 



eom^c 



(9) 



Here, we have approximated the synchrotron and self-Compton 
photons as isotropic in the jet frame, the average photon time in 
the source as 0jz/c, and the integration upper limit in equation ([8) 
accounts for the Thomson limit (assuring that e < 7) and 7 > 1. 
Equations ll8}j9]l can be solved iteratively, which yields all orders 
of Compton scattering. In the cases considered below, the first or- 
der dominates. For a synchrotron spectrum above the turnover en- 
ergy, e,, the integral can be performed analytically for the first-order 
Compton scattering. 



Jsc(e) ■ 



Ci{ajKB,,fc@,^R^ I B\~ 



967r2 



ficr 



min(e/7n, 1/e) 
In ° . (10) 
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dF _ 

AE ~ 


DJj(e)V 


0, 


1 


D^m^c'^ ' 


rj(l-/3jcos0' 



Then, we calculate the photon energy and the flux in the ob- 
server's frame. For the emission of the jet, we have 



(11) 



where E is the observed dimensional photon energy, Dj is the jet 
Doppler factor, and the fl ux transformation t o the observed frame is 
for a steady-state source (ISikora et al.lll997l) . The flux component 
from the counterjet should be then added, for which the Doppler 
factor, Dcj , corresponds to -ySj . Hereafter, we take into account the 
counterjet in numerical calculations, but since the jet emission for 
our adopted Cyg X-1 parameters (; = 27°, /Jj = 0.6) is higher than 
that of the counterjet by ~ (Dj/Ocj)- - 11, we neglect the coun- 
terjet spectral contribution in some flux estimates. Here, the form 
of dF/dE above assumes the energy unit in F and E is the same, 
whereas they are often assumed to be different (e.g, as erg and eV, 
respectively), which, however, can be easily accounted for. 

In equation Jill ). V is in the jet frame. A characteristic volume 
between z - z/2 and z + z/2 equals approximately 



V^n&j((R,)%. 



(12) 



This also equals the volume of a cylindrical region with the radius 
@j(Rg and height z = iRg. 

The optically-thin synchrotron emission from the considered 
jet and counterjet regions in the observer's frame is then given by 




10.0 



Figure 4. An example of the synchrotron jet spectrum (black solid curve) 
for p = 2.5, showing the spectral transition from the optically thick regime 
to the optically thin one. The normalization corresponds to the integral part 
of equation J2H . The dashed and dotted line shows extrapolation of the 
optically thick and optically thin regime, respectively. 



dFs _ CicrTK((R,)\&,B,,f l B 
dE ^ A'&nm^cD^ \fi,., 






(13) 



Around Ei, this spectrum joins the optically-thick emission. The 
optically-thick spectrum of Cyg X- 1 shown in Fig. [T] with a ^ Q 
has the flux of F, ^ 20 cm"^ s"' (== 13 mJy), which approximately 
agrees with the fits of Rl 1 . (We note that since the ~ 0. 1 eV region 
is strongly dominated by the emission of the star and its wind, see 
Fig-El the actual turnover energy may be different than ~0.1 eV.) 
We normalize then the optically thin synchrotron spectrum to that 
from the jet base, at zo, as dFs/dE(Et) = F,. Neglecting the coun- 
terjet, this yields the normalization of the electron distribution at 
the base, Kq, vs. Bq, where Bq is the value of B at zo. 



^ ^ 4S7TF,m,cD^ ^-ql /So 



F, 



(14) 



Then, a relationship between 60 and ^0 at the jet base 
is provided by the requirement that the synchrotron spectrum 
is self-absorbed below Ft = eiOjin^c^. The absorption optical 
depth in the observer's direction can be approximated as Ts - 
2q's(F)0jZo/-2)j sini, where 20jZo/sini is the path length going 
through the jet spine in the observer's frame, and Dj sin i is the 
sine in the jet frame. Taking into account equation 04), we solve 
Ts = 1 at F, for. 



i^Q- 



48FtC2mec\- ttD 



afCiQjSin// R,B^,\DjB„ 



So 



F, 



(15) 



This relationship is only weakly dependent of p through C2/C1. 
The inverse relation has Bq oc ^*. 

4.2 The partially self-absorbed jet 

We then calculate the synchrotron emission from the entire jet, for 
a conical jet without energy losses, as assumed by BK79, i.e., with 

K = Ko(aior\ B = Bo(a(or'- (i6) 

The source function for synchrotron radiation in the jet frame is, 

;s(6) C,a,cBli'B-^'^e''^~ 



5s(e) = 



Q's(e) 



24;r3C2 



(17) 
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In order to obtain the spectrum observed from the jet, we need to 
integrate the source functi on over the projected area of the jet. This 
implies, cf. equation (3) of lHeina <2006l1 



dE 



a' 



smj 

dzSs 



dz5< 



dx(l-e 



r 

j-0, 



d.Y(l-e-"s) + 



(18) 



where (for the jet) 
2as(f) 



Ts{€,Z,x) ■■ 



Dj sin ; 



(&izf 



ni/2 _ 2as(€)@jZ^ 
J ~ Djsin; 



(i-if^y\ (19) 



if/ = x/0jz, and 6 needs to be transfoiTned accordingly to the jet and 
and counterjet formula, respectively. As in Section \4~T\ we define 
the turnover energy by Ts (ft, Zo 1 0) = 1, i.e., at thebaseof thejetfor 
the path going through the spine of the jet. Then, 



TS 






^^ ' il-^)^ 



(1 



.^2)i 



(20) 



where ^ = (f/fo)(f/fi)- This allows us to express the flux observed 
from the jet as 



dFs 
dE 

X 



Q'fCi(Bc,<o^g)-0j sin/ /BjS, 



^1 



£d^f^£d^{l-exp[-r^(l-^^)^)] 



(21) 



An analogous term for the counterjet (for which sin(;r - /) = sin 
needs to be added, including the turnover frequency transfoiTned 
accordingly. In the optically-thick case, £ <c £,, and we can set the 
lower limit of the outer integration to 0. The resulting dimension- 
less double integral depends on p only, and we denote it as C3. Its 
values for p = 2, 3, 4 are ^ 3.61, 2.10, 1.61, respectively. Then, 
dFs IdE in the optically-thick case has o- = 0, 



F, 



afCiCAB^^oR^)-®, sin;- 10,8,, 



2A7t'C2m,cD- 



Bq 



1 + 






,(22) 



where Et is the jet turnover energy. We can see that this equation 
times n/ilCi) ~ 1 and neglecting the counterjet contribution is 
equivalent to equation dlSl l. This shows that the one-zone model 
provides a good approximation for the jet emission at Et (though, 
of course, the one-zone model does not yield the optically-thick, 
a = 0, spectrum). We can also substitute. 



p+i p+1 



ai sin ; 



Sc, 



(23) 



[which follows from rs(et, Zo,0) = 1] for £, in equation J21b to get 
the optically-thick flux expressed entirely through the intrinsic jet 
parameters. We note that, in a dependence opposite to the relativis- 
tic beaming along the jet axis of Ft oc 2) ''* '' ''* (for the jet), the 
optically-thick emission is also beamed away from the jet axis in 



the jet frame. Ft oc (sin ;)<''""'''''+'", and F, oc £)' 



,(p+2)/(,,+4) 
J 



/ sin /. 



Note that he transformed both the emission and absorption coefficients to 
the observer's frame by multiplying by D? {6^ in his notation) whereas the 
latter should be divided by Dj , as well as transfomiation of the frequency 
should be taken into account. Consequently, the powers of S in subsequent 
equations need be modified. Also, his emissivity is defined integrated over 
all directions and the magnetic pressure is assumed to be B-JlAn. 



Then, for E » Ft, ts(F) <c 1, which allows us to 
perform the double integration in equation (1211 analytically, as 
(F/Ft)"*''"'"^;r/(/7 - 1). We can then substitute Ft of equation ( |23l > 
to obtain the flux in the optically-thin case, 

dE lAnip - \)m,cD^ \b,,) \m,c^) \ J <^i j'^ ' 

We see that, for a given Kq, this equation times {p - l)/2 = a ~ \ 
equals equation ( I13t . which shows that the one-zone model pro- 
vides also a good approximation to the optically thin synchrotron 
emission. Note that the optically-thin emission is (in our approxi- 
mation to synchrotron) isotropic in the jet frame. 

Only in the transitional region we need to perform the integra- 
tion numerically. Fig.|4]shows the resulting spectrum (for either jet 
or counterjet) for p = 2.5. We see that the transition is rather grad- 
ual, not showing a sharp break, in contrast to the broken-power-law 
fit of Rl 1. We also note that the intersection of the optically thick 
and thin power laws, dFs/dF(thin) = Ft, occurs not at Ft (defined 
by Ts = 1 for the jet) but instead at a F°''% given by 



£. 



Clip - 1) 



(25) 



The ratio of equation ( I25t is == 0.75 for p = 2-5. We take into ac- 
count this correction to Ft in our models in Section l431 Neglecting 
the counterjet, the normalization of the electron distribution at the 
base becomes 



Ko 



lAjf-Fim^cD^ 



Bo 






O. 



(26) 



CiC,ajBl®]{(oR,f \Bc,y 

Note that this equals equation l ll4t multiplied by nlld. 

In order to determine the values of Bq and <fo separately, we 
can assume a degree of equipartition between the pressure or en- 
ergy densit y of the magnetic field and the relativistic electrons, 
as done by Ichatv, Dubus & Raichoon ( 1201 iHi for the black-hole 
binary XTE J1550-564. The usual quantity to describe it is the 
plasma parameter /?, which is the ratio of the particle pressure to 
that of the field. Since our calculations do not determine the pres- 
sure of ions, we define/? here for the relativistic electrons only, 

B^/Stt Bllin ' 

where we have assumed the magnetic pressure to be B'/Sn. A con- 
tribution from ions would change the actual value of/?. Here, u, is 
the energy density of the relativistic electrons, and / = ii^/Kin^c^, 
see equation ( 142 1 below. Given the dependence of K and fi on z of 
equation J16t . the value of j8 is constant along the jet. We can then 
substitute Kq of equation ( 126b in equation | |23I ). and solve it with 
equation ( 122b for 60 and <fo- Neglecting the counterjet, we find 



Bn - BcT 



-t 



O, 



lp+3 



io- 



E,R,^-^ 

3C2 
Cf sin i 



(lof sin ifCiCim^c^ f- 



C1C3 

2/J+13 



(37rC2)3(/?S„o-TO)20j/^, 



(28) 



p+6 p+1 



&>Bl 



D. 



(29) 



^ Note that our/? and 0j equal their 3f and hT^ , respectively. Both the right- 
hand side of their equation (A2) and the term mcc/3e should be multiplied 
by In. They neglect the relativistic effects and thus their counterjet flux 
equals that for the jet. 
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Then, we consider the self-Compton component. In the BK79 
model, the monochromatic optically-thin synchrotron flux per In ^ 
is oc (j'-(p-')/2 On the other hand, the self-Compton flux in this 
model decreases with height one power of i^ faster, oc ^-(p+i)/2 jj^js 
is due to the dependence of the Compton rate on K-, compared to 
oc K for synchrotron. Thus, this component is dominated by the 
base even more than the optically-thin synchrotron. Therefore, in- 
stead of solving the self-Compton emission from the entire jet, we 
consider here only Compton scattering of synchrotron photons pro- 
duced in the base region. 

The self-Compton component is most important at high en- 
ergies, where it dominates over the synchrotron one. At those en- 
ergies, its emission is dominated by the cooled part of the elec- 
tron distribution (|2), with the index p 4- 1 at 7 > yj,. Neglecting 
the high-energy cutoff and the counterjet emission (which we take 
into account in the numerical calculations in Section l43) . the self- 
Compton flux can be calculated using equations l |10J41U . 



p+4 



d£ 



X In min 



96mnscD^ 
D,Boyl ' EBorl 



(30) 



In the present formalism, we assume the turnover energy corre- 
sponds to the un-cooled electrons with the index p; then A^o is given 
by equation l |23t . However, it may also correspond to the cooled 
electrons, in which case equation ( I23t would need to be modified 
(by substituting p+ 1 foi p and Kojb for Kg). 

The main observational constraint on the amplitude of the self- 
Compton component comes from observations at ^ 0. 1 GeV. We 
use here the AGILE upper limit, Fy(Ey), as given in Section|2] Since 
the self-Compton flux decreases with both the increasing magnetic 
field and the increasing height along the jet, the constraint, 

dpQc 



dE 



(31) 



translates into a lower limit on Bo for given fo and a lower limit 
on <fo for given Bq. The values of 60 and fo are mutually related by 
equation l l22t . Using this, equation OU yields the limit for Bq, 



'Bo 

[ BcT 



Ftrlafsmi I £, 



AFyC2Cin 



m^c- 



m^c^ 



- 9 

In ' \ , (32) 
EyBorl 



independent of 0j and D. The corresponding limit on (q can be 
readily obtained using equations ( I30ti22t . Note that equations l l30t 
and l l32t depend on 71,, which itself depends on Bq, and the above 
limit depends logarithmically on Bq. Thus, these equations need to 
be solved iteratively. 



4.3 The jet parameters 

Given the knowledge of the electron distribution and the magnetic 
field, we can calculate the components of thejet+counterjet power, 

P. = yMe/3JC(^j0J^i^g)^ (33) 

Pp = 2nji^n,m,c'/3,T,(T, - l)(0jf/^g)^ (34) 

Pb = ^/S,c(r,@,(R,)\ (35) 

where P^, Pp and f g, is the kinetic power (i.e., the energy flux car- 
ried by the jet) in the electrons (including any positrons and taking 
into account the contribution from pressure), protons (assuming to 



be cold and not including the rest mass), magnetic field, respec- 
tively. Then, u^ and n^ are the energy density (including the rest 
mass) and the number density, respectively, of the electrons respon- 
sible for the observed emission, and zyp is the number of protons per 
emitting electron. Note that T/p can be < 1 if there are positrons in 
the jet, or > 1 if not all electrons in the jet are accelerated into 
the distribution (|2}. The powers of equations (I33ti35) can be cal- 
culated at any point of the jet, and their constancy with ( follows 
from equation ( |16t . 

In order to calculate ;/(., «(., and the synchrotron power, we 
need to know the full distribution of electrons. The electron adi- 
abatic and synchrotron (neglecting self-absorption) loss rates, and 
the Lorentz factor at which the two rates are equal are. 



7ad 



3z 



-7. 



4 o-T B^ , 

Ts = T —7 

3 m^c on 



7s 






(36) 



respectively. The factor of 2/3 in y.^^i accounts for the expansion 
being in two dimensions only and we have assumed that the jet is 
conical. AppendixlAlconsiders in more detail the dependence of 7s 
and the corresponding synchrotron photon energy on parameters of 
accreting sources with jets. We then define the Lorentz factor of 
electrons emitting at the local turnover energy, E,((), 



7i - 



Scr£,(<r) 



BDjin^c^ 



(37) 



Note that electrons with 7 < 7t, emitting synchrotron emission be- 
low Et, have a synchrotron loss rate much lower than 7s, and their 
adiabatic losses will dominate (neglecting Compton losses). Con- 
sequently, the cooling break energy below which adiabatic losses 
dominate will be approximately at. 



7b - max (7s, 7,). 



(38) 



Note that equation l[T6j and E,(() cc ^-' (BK79) imply that 75 cc (, 
7( oc i^". Here, we neglect this complication and assume a constant 
value of 7b corresponding to the jet base. The steady-state electron 
distribution above 7b has the p + I index. 

We assume then a high-energy cutoff is due to radiative loss 
rate becoming faster than the acceleration rate. Equating 7s to the 
electron acceleration rate occurring on a gyroperiod time scale, 7acc, 



ccB 



7iK 



2nm,c ' 



(39) 



where j/acc ^ 1 is a scaling factor, the maximum Lorentz 
factor and the maximum synchrotron energy are obtained as 
JGuilbert. Fabian. & Reeslll983l : lde Jager et all 19961) 



7; 



9B„7/„ 



97/a, 



' 50;/a, 



(40) 



S/ro-fB Snai 

As noted first bv lGuilbert et al.l ( Il983h . fm is independent of B. In 
equation l|2]l, we assume e-folding at 7m/2, which, given our delta 
function approximation to the synchrotron spectrum, equation lO, 
results in an exponential cutoff of the photon distribution at 6,^. In 
order to have e,^ - 2 (neglecting relativistic energy corrections), 
as fitted in the two-component model to Cyg X-1 (Section l3.2t . 
'7acc - 0.05 is required, which value we assume. 

Thus, we now adopt the steady-state electron distribution of 
equation l|2j, keeping in mind that the actual one may be more 
complicated. The main complication not taken into account in our 
presented formalism is the possible presence of a minimum accel- 
eration Lorentz factor > 1 (see Section l43t . With equation l|2j, we 
can calculate n^, u^ and Ps, the power of the synchrotron radiation 
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emitted in all directions by the jet and counteijet, 
K 



p-\ 



1 






Km^c^f, f : 



So 
Ps - — crTcA:oV/o 
3n 
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+ rb 



3-p 

Tb ■ 



„3-/' 






r 1 



2yb 



m 27b 



(41) 



(42) 



(43) 



erg s"'. The resulting GeV flux is > 3 orders of magnitude below 
the observational upper limit. The model yielding the GeV flux at 
the upper limit, shown in Fig.|5jb), has (o - 0.8 x 10^, So == 3 x 10^ 
G,/3 ^ 480, rs ^ 780, y, ^ 48, and („, ^ 2600. Also, Pe ^ 3 x lO^"* 



ergs 



2x10' 



ergs 



,xlO 



'3'ergs-',Ps == 1x10^'' 



where Vo is the volume, equation J12t . at fo, / is the dimensionless 
energy density, T is the incomplete gamma function, and we have 
assumed y^ » yb. In equation l l41t . we have also neglected the 
contribution of the high-energy cutoff, and in equation ( 143 K the 
emission below £,. We note that Pg has been integrated over the 
entire jet length [with the dependencies of equation (I16H. but it 
happens to be equal to the power emitted by the two base volumes 
in the one-zone approximation. 

Finally, we consider the constraint Pb = JJbP], where we use 
the observational estimate of Pj, and 77s < 1 from energy conserva- 
tion. Using equations ( 122b and ( |35l >, we obtain 



(0<^m 



/3jO-T 



T" n TlT^ m^ CT Di 



EtBa&- Rg 



(TjDj)? 



3F,C2 



fffCiCs sini 



(44) 



and (0 = imri's 



1/10 



erg s~ . The value o f Pp is only slightly above the range of Pj of 
iRussell et alj J2007I) . and it can be within it for a slightly higher 
value of Bo (yielding a lower GeV flux). We note that the weak 
magnetic field of this model, with /3 » 1 , is similar to that found 
in the y-ray emitting region of the jet in Cyg X-3 (Z12). In both 
models 1 and 2, the jet in Cyg X-I appears not dominated by e* 
pairs. 

We note that the acceleration rate may take place only above 
certain minimum Lorentz factor, yi » I. In Z12, yi ~ 10^^ was 
found necessary to explain the y-ray spectrum of Cyg X-3 from 
Fermi dFermi LAT Collaborationll2009l) . We have thus also consid- 
ered a model with yi = 10' . The model 3 accounts for the MeV tail, 
and it has p = 1.6. The spectrum has an o- == 0.5 power law above Et 
(due to electrons at < yi which are efficiently synchrotron-cooled), 
breaking at ~ 0.5 keV to a - 0.8. Given the relative similarity of 
the spectrum to that of our model 1, we do not show it. Models with 
yi » 1 and p > 1.6 are also possible. They do not explain the MeV 
tail, and are relatively similar to our model 2. 

Our models 1 and 2 have ys > yt, which appears compati- 
ble with observations of relatively hard optically-thin spectra abov e 



This limit is only weakly depending on p through £, JCorbel & Fendedl2002l : iGandhi et al.]|20Il(: IChatv et al.lboilh 



S tically -i 
mile 



4.4 Application to Cyg X-1 

Here, we consider the full-jet formalism of Sections 14.2114.31 In 
our model 1, we reproduce the MeV tail by the jet synchrotron 
emission. Given that the electrons responsible for X/y-ray emission 
are cooled, y > y^, (see Appendix |A), we have first considered a 
model with p = 1.2, at which the slope after the cooling brake is 
2.2 and the X-ray energy index is cc = 0.6. However, we have found 
that although the X-ray-emitting electrons do have y > yb, there is 
a hard [a = (p- l)/2 = 0.1] part of the spectrum above Et emitted 
by p = 1.2 electrons because yb » yt- This results in the I-MeV 
flux being above that observed. We have thus increased p to a value 
at which the model would match the I -MeV flux. Such a model has 
p = 1.35, and it is shown in Fig.|5ja)- 

We have first imposed equipartition between the electron and 
magnetic pressure, /? = I, see equations ( 12814291 ). This model yields 
the GeV flux somewhat above the observational upper limit, and a 
slight increase of 60, corresponding to/3 == 0.8, gives the GeV flux 
at that limit, see Fig.^a). The model has (0 ^ 860, Bo - 8.7 x lO"* 
G, f„ =: 2100, yb = ys = 80 > y, == 20. The low value of p 
of the model results in a rather slight break at Et, from a = to 
a = 0.175, which may be not compatible with the break found by 
Rll. The spectrum steepens to o- = 0.675 only around 1 eV 

The model has Pe == 8 x lO'"* erg s"', P^ ^ 2 x lO'"* erg s ', 
Pp == 777p X 10^'* erg s"' , Ps ^ I X 10"*^^ erg s"' . We note tha t P„/rj„ is 
much less than the estimated jet kinetic power, Pj, of Rus sell et al.l 
( 120071) ■ which requires the presence of j/p ^ 10^ times more protons 
than those coiTesponding to the emitting electrons. Increasing (0 
(allowed by the GeV upper limit) decreases P^i, Pp and /?, while it 
increases Pg and it does not affect Pg. 

In our model 2, we assume p = 2.5, motivated by p ^ 2.5 ±0.5 
found above £, in the black-hole binary GX 339-4 ( IGandhi et al.l 
I2OI ih . Imposing equipartition. B = I, we obtain <fo - 1.2x10^, Bo - 
I.lxlO'* G, similar to those found for model 1, and Pp - 377pXlO'^ 



However, if jet magnetic fields are strong enough (which, in case of 
Cyg X-1, is compatible with its GeV upper limit), ys < yt is pos- 
sible, which is indeed the case for our model 3. Then, the turnover 
region will correspond to cooled electrons. In this case, the elec- 
tron index used in Sections |4. 1114. 21 (which treatment is based on 
the synchrotron emission and absorption around £,) equals p + I, 
see equation ([2](, which will require certain modifications of our 
presented formalism. 



4.5 Irradiation of the jet 

The jet will be irradiated by the central X-ray source, including 
an optically-thick accretion disc, and by the star (see Section |2ll. 
The Doppler factor of the central source with respect to the jet is 
Dx = l/[rj(l -ySj)] - 2. The synchrotron losses will dominate over 
X-ray Compton losses (neglecting Klein-Nishina correction) for, 



B>Bx: 



(2Lx/cy''- 



ix 10- 



m 



G. 



(45) 



Our models have comparable values of Bo, and this effect should 
be taken into account in a more detailed treatment. Note that for 
the standard dependence of B oc z^' above the jet base, the ratio be- 
tween the X-ray Compton and synchrotron losses will be constant 
along the jet. 

The Compton-scattered emission will be mostly directed back 
towards the X-ray source due to anisotropy of relativistic Comp- 
ton scattering, and only a small fraction, a few per cent of the 
average Co mpton emission, will b e directed towards the observer 
at; ^ 27° dijhisellini et alj|l99ll : lOubus. Cerutti & Henril IJoIol. 
IZdziarski et al.l 120121) . The direct emission will contribute to the 
high-energy soft y-ray tail observed in Cyg X-1, and the one re- 
flected from the accretion disc will contribute to the observed 
Compton reflection spectral component. On the other hand, the 
same kind of emission from the counterjet will be beamed along 
the jet axis, but it will be blocked by the disc. Its absence may put 
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Figure 5. The data are the same as in Fig.[T] The red sohd and blue dashed curves show the model synchrotron and Compton components, respectively, (a) 
The model 1 for p = 1.35, accounting for the observed MeV tail, which corresponds to the approximately maximum jet emission allowed by the data, (b) The 
model 2 for p = 2.5, in which case the jet emission is well below the MeV tail, with it most likely being emitted by hybrid plasma in the accretion flow. See 
Section l4^ for details. 



some constraints on the outer radius of the disc (which may be rel- 
atively small in wind-fed accretion, which appears to take place in 
CygX-1). 

The magnetic field which energy density equals the stellar 
photon energy density at the height z along the jet is, 



B. 



(2L,/c) 



1/2 



D. 



1 



Fjd -piz/R) 



R = (z^+a^y'\ (46) 



At z <K a, S, =: 110 G, which implies the stellar irradiation at the 
jet base is completely negligible. We note, however, that for B oc 
z"' above the jet base, the Compton losses on the stellar emission 
may become dominant at large heights. For our models 1 and 2, 
the stellar Compton losses overcome the synchrotron losses at z ^ 
0.05a, 0.02a, and become 60 and 600, respectively, times stronger 
than the synchrotron losses at z ^ a. Thus, this appears to be a 
highly important effect, which should be taken into account for an 
outer jet of Cyg X- 1 . 



5 DISCUSSION AND CONCLUSIONS 

In Section[3l we have calculated the average hard-state spectra from 
INTEGRAL ISGRI and PICsIT detectors over the period of 2003- 
2010. They are found to agree well with those from CGRO (M02). 
This is consistent with the long-term X-ray and radio behaviour of 
Cyg X-1 in the hard s tate being rather constant since at least 1995 
IZd ziarskietal.l201ll) . and, in hard X-rays at > 20 keV, since 1991 
(Zdz iarski et al.l20 02il . 

We confirm the presence of a high-energy tail in the == 0.5-5 
MeV range. We find this tail to be very similar to that measured by 
CGRO, but much weaker th an that claimed by LI 1 . This has been 
independently confirmed bv ljourdain et al.l ( l2012h . and then by the 
revised measurement of Ph. Laurent (private communication). 

The origin of the tail is not clear. Although we find it is well 
fitted by hybrid Comptonization (presumably within a hot accretion 
flow), it can also be due to a separate spectral component, in partic- 
ular synchrotron jet emission provided it has a high-energy cutoff 
at ~ 1 MeV. If the measurement of the strong polarization above 0.4 
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MeV of Lll is confinned, the latter interpretation has to hold. We 
note, however, that the polarization fraction of 0.67 + 0.30 (Lll) 
is difficult to explain even by synchrotron models, as it requires 
extremely well ordered magnetic fields. For example, polarization 
fraction of synchrotron radiation in blazars never exceeds 50 per 
cent, it very rarely reaches 40 per cent , and the typical va lues are 
of the order of ~10 per cent (see, e.g.. |jorstad et al.ll2007l and ref- 
erences therein). A polarization fraction of only ^2A per cent was 
obse rved in the IR emi ssion of the black-hole binary XTE J 1550- 
564 JChatv et al.ll201lh , interpreted as optically-thin synchrotron 
jet emission. 

We note here that the electric vector position angle (hereafter 
EVPA) found by Lll is 140° ± 15°. They claim it is at least 100° 
away from the position angle of the radio jet in Cyg X-1, implying 
the polarization is perpendicular to the jet. However, as we noted in 
Secti on[2l the actual pos i tion angle of the jet in Cyg X-1 is -(17°- 
24°) l lStiriing et al.l200ll:lRushton et 31120111) . Then, given that the 
EVPA is invariant upon adding ±180° (since the electric field of 
a photon oscillates, changing its sign), the EVPA of Lll can also 
be written as -(25°-55°), and thus it is approximately consistent 
with the observed jet angle. Polarization parallel to the jet can be 
produced if the magnetic field is dominated by the toroidal com- 
ponent, and/or due to compression of chaotic magnetic field in the 
internal shocks formed with fronts oriented (quasi-)perpendicularly 
to the jet axis. 

In order to test whether the MeV tail may be the high-energy 
end of the optically-thin jet synchrotron emission, we have devel- 
oped a jet model allowing us to determine the parameters of its 
source. We noted that the synchrotron power-law spectrum begins 
at the turnover energy (Ef), which corresponds to energy at which 
the base of the jet becomes optically thin. This implies that most of 
the optically thin emission originates in the base, which we model 
as a single zone. Combining the expressions for the synchrotron 
emission and absorption at £,, we find that the magnetic field and 
the height of the jet base satisfy a relation, Bq « Zq- Physically, it 
follows from the requirement of obtaining a given synchrotron flux, 
which yields the total number of electrons for a given Bg, combined 
with the requirement of self- absorption at Et. The higher up the 
emission originates, the lower the electron density (for a conical 
jet) and thus the lower self-absorption optical depth for a given Bq, 
which has to be compensated by an increase of Bq. 

We have then developed a formalism for the synchrotron emis- 
sion of the full jet (reformulating the model of BK79) at any energy. 
The emission is found by analytical formulae in the optically thick 
and thin regimes. We have confirmed the validity of the one-zone 
model for the optically thin emission. Our results give also detailed 
spectra in the transition region, near Et, as an integral. The spectra 
have a smooth transition between the two regimes, whereas Rl 1 
fitted the transition region as a broken power law. This may affect 
the accuracy of their values of £, and E(Et). 

We then consider synchrotron self-Compton emission. For a 
given synchrotron flux at £(, a monochromatic flux from this pro- 
cess is oc S(, . We use the observational upper flux limit at 
0.1-3 GeV from AGILE. This gives a lower limit on So, and using 
the found relationship between zo and Bq, also on zq. Next, noting 
that the jet Poynting flux c annot exceed the total jet kinetic powe r 
estimated observationally JGallo et al.ll2005l : iRusseU et al.ll2007h . 
we obtain upper limits on Bg and zq. We also calculate Bq and zo 
from the assumption of equipartition between the magnetic field 
and relativistic electrons. 

The above results are applied to Cyg X-1, using the values 
of El and F(Et) of Rl 1. This yields the location of the jet base at 



~ 10'i?g, and Bo ~ 10"* G, relatively weakly dependent on the as- 
sumed value of p. We also find that in order for the observed MeV 
tail to be due to the jet synchrotron emission, p ^ 1.3-1.6, which 
are relatively low values for acceleration processes. We note that 
such acceleration indices are harder than that assu med before in 
X-ray jet models of black-hole binaries, p - 2-2.3 |Markofr et al.l 
'2003'; 'Heinz & Sunvaev"2003'; 'Merloni, Heinz & di Matted 1 20031 : 
Falcke, Kording & Markoff 2004; Heinz 2004). Low values of the 
acceleration index appear also to be in conflict with observations 
of opt ically-thin synchrotron radio spectra, e.g.. lMiller- Jones et al.l 
( I2OO4I) . On the other hand, values of p > 2 are consistent with 
standard acceleration models and observations of optical ly thin 
synchrotron spectra (e.g., Bednarz & Ostrowski 1998; Kir k et al.l 
I2OOOI) . Still, our present understanding of electron acceleration ap- 
pears insufficient to rule out models with p <2. 

We find that models with p > 2 can have a much higher kinetic 
jet power, due to a higher number of ions associated with the emit- 
ting electrons with a steep spectrum, as well as due to the magnetic 
field allowed then to be low by the GeV upper limit. Our model 2 
has t he kinetic power ab out equal to that determined observation- 
ally JRussell et alj|2007r) . This class of models, however, have the 
synchrotron MeV fluxes much below those of the MeV tail, and 
thus the tail cannot be due to that process. This, and the spectral 
proximity of the tail to the high energy break of the hot disc emis- 
sion may favour production of the tail in a hot accretion disc, as 
predi cted b y the hybrid Comptonization model l| P outanen & Coppil 
i 19981 ; M02; |Poutanen& Vunnll2009l ; lMalzac &Belmontli2009h . 

These values of Bo and fo are comparable to those recently ob- 
tained for XTE J1550-564 assuming a magne tic-electron equipar- 
tition in a one-zone model JChatv et all 2011*) . and for GX 339-4 
JGandhi et al.ll201in . In these objects, the optically- thin emission 
above the turnover frequency shows p ^ 2-3 JCorbel & Fenden 
I2OO2I ; iGandhi et al.ll201 ll ; IChatv et aljboi ll) . Given the necessary 
cooling break (if those indices correspond to the un-cooled part of 
the electron distribution) corresponding to an energy between the 
IR and X-rays, the implied jet contribution to the X-rays is low. 
Interestingly, GX 339-4 al so shows a hard X-ray tail on top of 
thermal Compton spectrum (IWardziiiski et al.l2002l ; lDroulans et al.l 
I2OI0I) similar to that in Cyg X- 1 , which tail then has to have the ori- 
gin different than synchrotron emission of the jet. 

We note that the location of the jet base in Cyg X-1 at zo ~ 
10''i?g ^ 2 X 10' cm is in agreement with the observed large orbital 
modulation of the radio emission (which is due to orbital-phase 
dependent free-free absorption by the stellar wind from the donor), 
with, e.g., the depth of ^ 30 per cent at 15 GHz, which requires the 
bulk of the radio flux to be emitted at a distance comparable to the 
orbital separatio n, z ~ a ^ 3 x 10'^ cm JSzostek & Zdziarskil2007l ; 
IZdziarskill2012h . If the turnover frequency (emitted by the base) is 
Vt ~ (2-3) X 10" Hz (Rl 1), yt/(15 GHz) ^ a/zo, in agreement with 
the prediction of the z »: v"' dependence of BK79. 

On the other hand, both our estimate of the jet-base distance 
and the strong orbital modulation of the radio emission are not con- 
sistent with the interpretation in terms of the one-component BK79 
model of ~50 per cent of the 8.4 GHz e mission observed to be re- 
solved at the scale of ~ lO'^'/sini cm JStirling et al.ll2001ii) . This 
interpretation yields the l ocation of th e t = I region for 8.4 GHz 
emission at z ~ 10'"* cm JHeinzluOOfil) , which implies the jet base 
at Zo ~ 3 X 10'" cm ~ W*Rg, an order of magnitude above our es- 
timates and implying virtually no orbital modulation at 8-15 GHz. 
A likely solution of this problem is an occurrence of a secondary 
dissipation event at a large radius, see a discussion in IZdziarskil 
(2012) . Then the inner jet may still follow the BK79 model while 
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the resolved emission is due to secondary dissipation. This also 
solves the problem of the Cyg X-1 jet power being much lower 
than the power inferred from the ring n ebula, with the jet p ower 
for our model 2 about equal to that of Russell et al. (2007), and 
much higher than t he theoretica l estimate of Heinz (2006) . (We 
note, however, that iHeina l200q did not include in his jet power 
estimate the power of the ion bulk motion, which we find to be 
dominant, see Section 1431 ) We note that two dissipation regions 
at very different scales, one emitting IR to y-rays close to the com- 
pact object, and one emitting radio far away, are present in Cyg X-3 
(Z12). 

Finally, we point out that electrons in the jet base are subject to 
substantial irradiation by the central X-ray source. For z ~ lO^Rg, 
its energy losses will be dominated by Compton up-scattering of 
X-rays. We also point out that although irradiation by the stellar 
photons is negligible at the jet base, it appears to dominate the syn- 
chrotron losses at distances along the jet comparable and larger than 
the orbital separation. 
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APPENDIX A: ELECTRON ENERGY LOSSES AND THE 
JET MAGNETIC FIELD 

In this work, we have found that the relativistic electrons responsi- 
ble for optically-thin emission at £ » £, are efficiently cooled for 
models satisfying the observational constraints. Here, we consider 
the effect of synchrotron cooling more generally. The importance 
of the break due to the radiative losses for t he electron and pho- 
ton distributions in X-ray jet models (e.g., of lMarkoffet al. 20031: 
iFalcke et alj |2004) was earlier poin ted ou t by iHeinz & SunvaevI 
1 20031) , and then studied in detail bv lHeina (2004). However, they 
have not calculated the possible values of this break as a function 
of the possible magnetic field in the jet. 

The Lorentz factor at which the synchrotron losses equal the 
adiabatic losses, ys, is given b y equation (I36t . Sinc e B oc M~'^^ in 
most equipartition models JHeinz & Sunvaevll2003l : iMerloni et al.l 
|2003J), and z oc M f or a given (, ys in those models is independent 
of M. We note that y/yad is both the time scale for the adiabatic 
losses and the characteristic time scale a flow element of the jet 
spends at ~ z- Thus, for y > ys, the loss process has enough time 
to significantly steepen the electron distribution. If the electrons 
are reaccelerated locally, their steady state distribution steepens by 
unity with respect to the accelerated distribution. If the electrons are 
only advected from lower heights, a high-energy cutoif develops 
above ys. 

The characteristic observed photon energy corresponding to 
ys, see equation l |36l ) (assuming ys > yt), emitted at a height of ( is 



o-jr"^GM(^ 



(A5) 



The scale factor, tjb, implied by equations l lA2t and l lASt is inde- 
pendent of ^, 



riB- 



o 1/2 ■ 

2r pj 



(A6) 



In order to obtain the break energy for this model, this rjg should 
be inserted into equation ([A3}. Both cases rule out the presence of 
an un-cooled jet X-ray spectrum, unless the magnetic flux is ex- 
tremely small compared to the jet flux. In that case, however, the 
self-Compton emission would yield veiy strong y-ray fluxes, in vi- 
olation of, e.g., the GeV upper limit for Cyg X-1 (see Section|4j. 



m,c- 






(Al) 



We scale the magnetic field with respect to the jet power [as in 
equation ll44tl. 



-l3,c{T,@,(R.f = nsPy 



(A2) 



This yields the cooling-break energy in the observer's frame of. 



E, _ 3'/^D,/^^W^ 

1.05 X io-"~DjpJ'^r^®^i 

~ 3/2 3/ 

< Pi 



3/2 



GMq) 



1/2 



(A3) 
(A4) 



where m = M/Mq, ao = r^/o'^ is the Bohr radius, and r^ is the 
classical electron radius. We see that to achieve an un-cooled power 



law extending to hard X-rays, E\y 



in a jet with a moderate 



or high power requires a very low value of tjb- 

We can also relate the magnetic fi eld in the jet to that in th e 
hot accretion flow, as in the models of iHeinz & SunvaevI ( l2003h . 
iMerloni et aT] ( l2003h and lHeinzl ( l2004h . The hot flow is, most likely, 
described by some variant of th e advection-dominated accretion 
flow (ADAF.'Naravan & Yill995[) . That model is quite close to that 
of spherical accretion, and, for the sake of simplicity, we will relate 
the pressure in the latter model, Pq = m^nv^ (where n is the pro- 
ton density and I'lf is the free-fall velocity), to the pressure of the 
magnetic field, S/Stt = bPo, where & is a coefficient that includes a 
hot-flow magnetization, a departure of the pressure in the hot flow 
from that of spherical accretion, and an overall scaling between the 
hot flow pressure and the jet pressure. The jet is launched from 
within a hot-flow radius, rjiJg. We then assume conservation of the 
magnetic flux, S oc z"'. This yields. 
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